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The noise-induced pattern formation in a population dynamical model 
of three interacting species in the coexistence regime is investigated. A 
coupled map lattice of Lotka-Volterra equations in the presence of multi- 
plicative noise is used to analyze the spatiotemporal evolution. The spatial 
correlation of the species concentration as a function of time and of the 
noise intensity is investigated. A nonmonotonic behavior of the area of the 
patterns as a function of both noise intensity and evolution time is found. 

PACS numbers: 05.40.-a, 87.23.Cc, 89.75.Kd, 87.23.-n 

1. Introduction 

The study of models of population dynamics is one of the topics that in 
the last few years has revealed large interest ranging from chaos to spatial 
organization, covering several branches of physics and different disciplines 
such as biology, theoretical ecology, oceanography and medicine. In par- 
ticular the addition of noise in such mathematical models can be useful to 
describe the phenomenology in a realistic and relatively simple form. The in- 
troduction of noise in real systems gives rise to non trivial effects, modifying 
sometimes in an unexpected way the deterministic dynamics. Examples of 
noise induced phenomena are stochastic resonance 012]) noise-enhanced sta- 
bility 3 , temporal oscillations and noise-induced pattern formation [IJ [31121 • 
Biological complex systems can be modelled as open systems in which in- 
teractions between the components are nonlinear and a noisy interaction 
with the environment is present |71 . The dynamics of interacting biological 
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species can be successfully described by means of a set of Lotka-Volterra 
equations with the addition of a multiplicative noise and a diffusive term 
to take into account spatial extension of the ecosystem. Recently it has 
been found that nonlinear interaction and the presence of multiplicative 
noise can give rise to pattern formation in population dynamics of spa- 
tially extended systems |H1 IHl- The model studied in this work concerns 
the ecosystem composed of three interacting species: two competing preys 
and one predator. The area of the patterns of maximum density has been 
quantitatively evaluated, finding a nonmonotonic behavior as a function of 
time and as a function of the noise intensity. The site correlation between 
the species concentration as a function of time and of the noise intensity is 
also investigated. 

2. The model 

The starting point of our study is the simple discrete set of general- 
ized Lotka-Volterra equations describing a population of two preys and one 
predator. The dynamics of our spatially distributed system is therefore 
described by the following model of coupled map lattice (CML) JHI 

x^+i = ^^xl^{l - vxl^ - fTvl^ - azl^) + V^^,^", + ^E(^p " 

p 

= - ^Vi, - - + ^,^^^3 + ^ E(yp - 

p 

= ^.zI:,[-l+7«,+<,)]+^/9<,^",+^E(^;-^I:.), (2.1) 

p 

where x"^, y"^- and zf^^ are respectively the densities of preys x, y and 
the predator z in the site at the time step n. Here a and 7 are the 
interaction parameters between preys and predator, X, Y and Z are the 
white Gaussian noise variables with 

{Xit)) = {Y(t)) = {Z{t)) = 0, (2.2) 

{X{t)X{t + r)) = {Y{t)Y{t + r)) = {Z{t)Z{t + r)) = 6{t), (2.3) 

{Xit)Y{t')) = {X{t)Z{t')) = {Y{t)Z{t')) =0 V (2.4) 

q is the noise intensity, D is the diffusion coefficient, ^ and are scale 
factors. indicates the sum over the four nearest neighbors in the map 
lattice. The multiplicative noise in the above equations models the inter- 
action between the species and the environment. The boundary conditions 
have been established in such a way that no interaction is present out of 
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lattice. This means that for the four corner sites we have only two interac- 
tions and for the other 4x98 line-confined sites the number of interactions 
is three. 

We analyze the transient dynamics of the system with a time varying 
interaction parameter P{t) between the two preys 

(3{t) = l + e + r]cos{ujt), (2.5) 

due to the environment temperature. Here r] = 0.2, u = ttIO"^ and e = 
—0.1. The interaction parameter P{t) oscillates around the critical value 
/3c = 1 in such a way that the dynamical regime of Lotka-Volterra model for 
two competing species changes from coexistence of the two preys (/? < 1) 
to exclusion of one of them (/?>!) The parameters used in our 

simulations are ^ = 2; = 1; q = 0.03; Hz = 0.02; 7 = 205 and D = 0.1. 
The noise intensity q varies between 10~^^ and 10~^. With this choice of 
parameters the interspecies competition among the two prey populations 
is stronger compared to the intraspecies competition (preys-predator), and, 
therefore, both prey populations can stably coexist in the presence of the 
predator pTT] . 

By considering the infiuence of noise and spatial diffusion, we analyze the 
noise-induced pattern formation. Specifically the time evolution of the area 
of the patterns and the correlation r over the grid between two species, as a 
function of the noise intensity, are analyzed. The quantitative calculations 
of the site correlation between a couple of species in the lattice have been 
done using the following formula 

EM(<,-^")(fci:,-fc") 
[ES«,-«)")2ES(fei:,-^")2 

where N is the number of sites in the grid, the symbols w'^,k'^ represent 
one of the three species x,y,z, and u;",A;" represent the mean values of 
the concentration of the species in all the lattice at the step n. Prom the 
definition 1)2. 6|1 it follows 

- 1 < r" < 1. (2.7) 

We calculate the correlation between the two preys and between preys and 
predator. 

3. Results and Comments 

The bidimensional spatial grid considered is composed by 100x100 sites 
in (x, y) plane. The calculations have been done for various noise intensities 
and at different steps of the iteration process. To quantify our analysis we 
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consider only the maximum patterns, defined as the ensemble of adjoining 
sites in the lattice for which the density of the species belongs to the interval 
[3/4 max, max], where max is the absolute maximum of density in the 
specific grid. For each spatial distribution, in a temporal step and for a given 
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Fig. 1. Spatial configuration of the species obtained after long time evolution in 
noiseless dynamics and in the exclusion regime between preys. The parameter set is: 
e = 0.1, /3 = 1.1,7? = 0,q = 0,D = 0.1, = 2,i/ = l,a = 0.03, /i^ = 0.02,7 = 205. 
The initial conditions are random with a Gaussian distribution, with mean values 
x(0) = y{0) — z{0) — 0.25 and variance (Tq — 0.1. Here ri2,ri3,r23 and ri23 are 
respectively the site correlations between: (i) preys, (ii) prey 1 and predator, (iii) 
prey 2 and predator, and (iv) predator and both preys. 

noise intensity value, the following quantities have been evaluated referring 
to the maximum pattern (MP): mean area of the various MPs found in the 
lattice and spatial correlation r between two preys, and between preys and 
predator. 

3.1. Deterministic Analysis 

In the absence of noise and with constant interaction parameter /? we 
obtain: (i) for e < (/? < 1) a coexistence regime of the two preys character- 
ized in the lattice by a strong correlation between them with the predator 
lightly anti-correlated with the two preys; (ii) for e > (/3 > 1) wide 
exclusion zones in the lattice (see FigQ, characterized by a strong anti- 
correlation between preys. In Fig^we report the spatial configuration of 
the species after long time evolution. We chose Gaussian initial distribution 
with mean value Xq = Vo = Zq = 0.25 and variance do = 0.1. 

By considering the periodic variation of the interaction parameter P{t), 
we obtain for e = 0, after a transient anti-correlated behavior between preys, 
a coexistence regime with strong correlation between preys that evolves to- 
wards an homogeneous spatial distribution of all three species. For e > 
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Fig. 2. Site correlation parameter r in noiseless dynamics as a function of time for 
different values of the parameter e : —0.1, 0., -f 0.1. Here 77 — 0.2. The values of the 
other parameters are the same used for Fig^ 
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Fig. 3. Site correlation parameter as a function of time. Here e = — 0.1,?7 — 0.2 
and the noise intensity is q — 10"^. The values of the other parameters are the 
same used for Fig^ The initial spatial distribution is homogeneous and equal for 
all species, i. e. x-™* = j/j"'* = 2:*™* = 0.25 for ah sites (i, j). 
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we find an oscillating behavior of the site correlation coefficient from co- 
existence regime between preys, corresponding to strong correlation, to an 
exclusion regime, corresponding to strong anticorrelation. This last be- 
havior is prevalent. The oscillating frequency coincides with that of the 
/3-parameter. When e < 0, the two preys, after an initial transient, remain 
strongly correlated for all the time, in spite of the fact that the parameter 
f3{t) takes values greater than 1 during the periodical evolution. This situa- 
tion corresponds to a coexistence regime between preys. In FigElwe report 
the behavior of the site correlation parameter r as a function of time for 
three values of the parameter e = —0.1, 0, 0.1. 



a) All the Species . b) Predator Only 
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Fig. 4. Semi-Log plot of the mean area of the maximmii patterns for all species as 
a function of noise intensity, at iteration step 1400. Here circles and triangles are 
related to preys, squares to predator. The values of the other parameters and the 
homogeneous initial distribution are the same used in Fig|2| 



3.2. Noise effects 

To analyze the effect of the noise we focus on the interesting dynamical 
regime characterized, in absence of noise, by coexistence between preys in all 
the period of /?, i. e. with e < 0. The noise triggers the oscillating behavior 
of the site correlation r giving rise to periodical alternation of coexistence 
and exclusion regime. Even a very small amount of noise is able to destroy 
the coexistence regime periodically as we can see from FigEl where the corre- 
lation parameter r as a function of time is reported. Noise is also responsible 
for a nonmonotonic behavior of the area of spatial patterns, which repeats 
periodically in time. In Fig0]we report a nonmonotonic behavior of the area 
of the maximum pattern as a function of noise intensity. A maximum of the 
area of maximum patterns is visible for the preys at g = 10~^ and for the 
predator at = 10~^. The same behavior is present in the following time 
steps within the first period of the interaction parameter: 600, 800, 1200, 
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Fig. 5. Spatial Pattern formation for preys and predator, at time iteration 1400 and 
for the following values of the noise intensity: q = 10~^^, 10~^, 10~*, 10"'', 10~^. 
The values of the other parameters and the homogeneous initial distribution are 
the same used in FigOl The parameters ?'i2, ris, r23, ri23 have the same meaning 
of FigUl 



1400. But at time steps 600, 800 the preys are highly anticorrelated with 
site correlation parameter ri2 = —1, while at time steps 1200, 1400 the preys 
are highly correlated with ri2 = 1. The formation of spatial patterns appear 
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only when the preys are highly correlated, while large patches with cluster- 
ization of preys appear when they are anticorrelated. This means that the 
coexistence regime between preys corresponds to the appearance of spatial 
patterns, while the exclusion regime corresponds to clusterization of preys. 
The noise-induced pattern formation relative to the iteration 1400 is visible 
in Fig|SJ where we report five patterns of the preys and the predator for the 
following values of noise intensity: q = 10'^^ , , W'^ , 10'^ , 10''^ . The 
initial spatial distribution is homogeneous and equal for all species, that is 
^vnit _ yinit _ ^tnit _ q 25 fgj- a,ll sitcs A spatial structure emerges 

with increasing noise intensity. At very low noise intensity, with respect 
to the value of the diffusion coefficient D, the spatial distribution appears 
almost homogeneous without any particular structure (see FigEK)- Increas- 
ing the multiplicative noise intensity, the symmetric dynamical evolution of 
the two preys in each site of the lattice is destroyed, so oscillations in pop- 
ulation density produce an exclusion regime of one of two preys [2]. This 
time evolution scenario corresponds to the appearance of spatial patterns 
due to different spatial density in each site. This spatial pattern disappears 
for sufficiently large noise intensity, producing a random spatial unhomo- 
geneity (see Fig|S^). As a final investigation we analyze the behavior of 
the area of the patterns as a function of time. We observe a nonmonotonic 
behavior of the area of MPs as a function of time for all values of the noise 
intensity investigated. Particularly for noise intensity values greater than 
q = 10"'' this nonmonotonic behavior becomes periodical in time with the 
same period of P{t), as shown in Fig |^ for q = 10~^. We note that this 
nonmonotonic behavior doesn't mean that a spatial pattern appears, like 
that of Fig|Sj3, but that a big clusterization of preys may occur, as shown 
in Fig[7J In this figure in fact we report the spatial configuration obtained 
for a noise intensity value q = 10~^, which corresponds to the maximum of 
the behavior of the mean area of MPs of preys as a function of time shown 
in FigE^. We see that this maximum corresponds to large patches of preys 
in the lattice investigated. The various quantities, such as pattern area and 
correlation parameter, have been averaged over 200 realizations, obtaining 
the mean values shown in the figures (3), (4) and (6). The effects induced 
by the interaction between the species and the environment, modelled by 
the multiplicative noise, can be summarized as: (i) to break the symme- 
try of the coexistence regime between the preys, producing an alternation 
with the exclusion regime; (ii) to trigger the oscillating behavior of the site 
correlation coefficient; (iii) to produce a nonmonotonic behavior of the pat- 
tern function of the noise intensity with an appearance of spatial 
patterns. 
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Fig. 6. Mean area of Maximum pattern of tlie tliree species and relative sites cor- 
relations between preys and between preys and predator as a function of time and 
q = 10"''. (a): black and white circles are related to preys, triangles to predator; 
(b) site correlation ri2 (black circles), ri3 (white circles), and ri23 (triangles). The 
values of the other parameters and the homogeneous initial distribution are the 
same used in FigEl 




Fig. 7. Spatial configuration of one prey obtained for q = 10~^, at iteration step 
800, corresponding to the maximum value of the mean area for preys reported in 
FigEl The values of the other parameters and the homogeneous initial distribution 
are the same used in Fig|31 



4. Conclusions 

The noise-induced pattern formation in a lattice of three interacting 
species, described by Lotka-Volterra generalized equations, has been inves- 
tigated. We find nonmonotonic behavior of the mean area of the maximum 
patterns as a function of noise intensity. The same behavior we find for the 
area of the patterns as a function of evolution time. The noise changes the 
dynamical regime of the species, breaking the symmetry of the coexistence 
regime. Besides the noise produces spatial patterns and temporal oscilla- 
tions of the site correlation parameter defined on the lattice. We finally 
note that our simple model of an ecosystem of interacting species could be 
useful to interpret the experimental data of population dynamics strongly 
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affected by the noise [T^ . 
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